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Abstract 

Besides their appearance at short distances > 1/M\y, local dimension-eight opera- 
tors also contribute to kaon matrix elements at long distances of order > l//i ope , where 
/K ope is the scale controlling the Operator Product Expansion in pure QCD, without weak 
interactions. This comes about in the matching condition between the effective quark 
Lagrangian and the Chiral Lagrangian of mesons. Working in dimensional regulariza- 
tion and in a framework where these effects can be systematically studied, we calculate 
the correction from these long-distance dimension-eight operators to the renormaliza- 
tion group invariant Bk factor of K° — K° mixing, to next-to-leading order in the 1 /N c 
expansion and in the chiral limit. The correction is controlled by the matrix element 
(0\slG ^ul^ 1 dL\K°) , is small, and lowers Bk- 



1 Introduction 



In order to study the physics at the scale of the K meson one can integrate out the fields 
describing all the other heavier particles. In the Standard Model, this integration generates 
a local AS = 2 operator of the formpQ 

«a/=2 = ^T^i + A?F 2 + 2X c \tF 3 ] C as=2 (h)Qas=2(x) (1) 
where Cas=2 is the Wilson coefficient given (in the large- N c limit) by 



Ca5=2(m) 



a s (n) / 1433 ac 
7T 1 1936 8 



9 

HJVc 



(2) 



The Wilson coefficient encodes the ultraviolet physics which has been integrated out and is 
computed in perturbation theory. Consequently, it depends on the renormalization scale (/x) 
and scheme (ac). In Eq.© k = or — 4 depending on whether one has used naive dimensional 
regularization or the 't Hooft-Veltman scheme (respectively) in the calculation. The operator 
Qas=2{x) in Eq. © reads 

Qas=2(x) = s L {x)^^d L {x)s L {x)^d L {x) , = 2 75 il> ■ (3) 

The Fi are functions of the heavy particles which have been integrated out and their expres- 
sion can be found in pQ. 

The invariant bag factor Bk is then defined in terms of this effective weak hamiltonian 

as 

(k \C as =2^)Qas=2(0)\K ) = ^F 2 K M 2 K B K (4) 

and, by construction, it is renormalization scheme and scale independent. 

Because of confinement and chiral symmetry breaking, there is another effective de- 
scription which involves directly the lowest-lying mesons as degrees of freedom in a chiral 
Lagrangian organized in powers of momenta and meson masses. In the case of AS = 2 
interactions, the lowest-order operator in this chiral Lagrangian is given by 

^A5=2 = ^A5=2(/i)Tr(A3 2 C/ t ^[/A 32 ?7 t r / ,C/) , (5) 

where (A32)y = <5?,3<% is a (spurion) matrix in flavor space, F ~ 0.087 GeV is the pion 
decay constant in the chiral limit and U is a 3 x 3 unitary matrix collecting the Goldstone 
boson excitations and transforming as U — > RU L* under a combined unitary flavor rotation 
(R,L) of SU(3)r x SU(3)l- In © <?AS=2(/-0 is a low energy coupling to be determined, as 
usual, by a matching condition to the Lagrangian (|T|). This condition is a crucial ingredient 
in the construction of the chiral Lagrangian ©and secures that the ultraviolet behavior of 
© be the same as that of © so that physics is independent of conventions such as the 
renormalization scale fx, the renormalization scheme (i.e. NDR vs HV), etc... 

Although matching conditions are commonplace in the construction of effective Lagrangians 
such as © in a perturbative context, it is only recently that their crucial role for the con- 
struction of the Lagrangian © has been appreciated also in the nonperturbative context of 
analytic calculations of electroweak matrix elements of light mesons [3]. 
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In the particular case of Bk and 3as=2(a*), a convenient Green's function to do the 
matching with is|2j 

W^ RI f(q, I) = lim i 3 J d 4 x d 4 y d A z e iqx e a ^ < 0\T{Lf(x)Ri s (y)Ll d (0)Rf{z)}\0 > (6) 



where 

Lf{x) = S^J-^d^x) Rf(x) = d(x) 7M i±^ a (x) (7) 

are QCD chiral currents. It is a general property that chiral coupling constants such as g&s=2 
in Eq. (JSJ) can always be defined by means of QCD Green's functions -in the forward limit- 
which are chiral order parameters, and hence vanish to all orders in perturbation theory in 
the chiral limit. 

Performing a standard analysis of Ward identities leads to the following relation 



dn q gr lUfctonzed = [ - g aP ) w LRLR {Q 2 ) , Q 2 = - q 2 , (8) 



where J d£l q stands for the average over the momentum q in D dimensions, namely 

dn q q,q u f(q 2 ) = ^ f(q 2 ) , (9) 
for a given function f(q 2 )- Using then the definitions 

z = -^- , W LRLR (z) = -^I-W(z) , (10) 

Vhad z ^had 

one can express gAS=2 in Eq. (JSJ) asjH 



32^ 2 Wlad) r(2-f)7o 



9AS=2^,e) = l-^f^[-^) I dzz~W{z). (11) 



This equation expresses the matching condition between the chiral meson Lagrangian (JSJ) and 
the quark Lagrangian We emphasize that the scale ^ 2 iad in Eq. (|TT|) is used solely for a 
rescaling of the momentum Q 2 and, in principle, is totally arbitrary. In practice, however, 
there is always a residual dependence on /U^ ad in the final result since calculations are done 
to a finite order in, e.g., a s . Therefore this scale fJ>i ac i has to be numerically large enough so 
as to make meaningful the truncated series in a s . Physically, one can think of //| a ^ as the 
scale at which resonances are integrated out and the Lagrangian © sets in. 

^From Eq. (Jll[) it is clear that g&s=2{lj) requires the knowledge of the Green's function 
W(z) over the full range of momenta and regretfully, we hasten to say, this information is 
not available. What is available, however, is its high- and low-z behavior since they are given 
by the Operator Product and Chiral expansions, respectively. Also, it is known that the 
analytic structure of W(z) simplifies notably in the large- N c limit becoming a meromorphic 
function, i.e. consisting of only poles, with no cuts. However, lacking the solution to QCD at 
large N c , the problem is still too difficult to tackle even in this limit: the masses and residues 
of the infinite number of resonances contributing to W (z) are not known. Obviously, a 
further approximation is necessary and this has been termed "The Hadronic Approximation 
to large- iV c QCD" (HA). 
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In order to motivate this approximation, it is important to realize that it is not the point- 
wise behavior of W(z) that is necessary, but only its integral. Moreover, W(,z) is defined fully 
in the euclidean regime. Therefore, it is reasonable to expect that a smooth interpolation 
between the low-z and the high- z regimes of W(z) may do a good job in approximating the 
integral. Consequently we shall construct the HA as a rational approximation 1 to large- N c 
QCD by keeping only a finite number of zero-width resonances whose masses and residues 
are obtained from matching to the first few coefficients of the Operator Product and Chiral 
expansions. In the case of the coupling g^s=2 in Eq. these expansions correspond to 

the high- and low- z behavior of the function W(z), respectively. 

The rational approximant so constructed, Wha( z ), is therefore an interpolating function 
between the high-z and low-z tails of W(z) which it matches by construction. The more 
terms one knows in the Operator Product and Chiral expansions, the more resonances one 
can determine from the matching to construct the interpolator Wha(z)- Therefore the ap- 
proximation is, in principle, improvable. In practice, however, often only one or two terms in 
each expansion are known, so that one actually has to borrow the masses of the resonances 
from the Particle Data tables and leave only the residues to be determined by these matching 
conditions 2 . There is a minimal requirement, though, which is that Wha{ z ) has to reproduce 
at least the leading non-vanishing term in the OPE. This is indispensable for accomplishing 
cancelation of scale and renormalization scheme dependence with Wilson coefficients, such 
as Ca5=2(m). This was done in the calculation of Bk hi Ref. 2]. Several other observables 
have also been computed within this framework[3] 0. 

We may now take the example of the matching condition (jllj) to discuss the appearance 
of local dimension-eight operators at relatively long distances. Indeed, the large-z fall-off of 
the function W(z) in this matching condition between the chiral meson Lagrangian and the 
quark Lagrangian is given by the OPE of the Green's function Dimension-eight operators 
(and higher) contribute to this OPE and thus they appear in kaon matrix elements at rather 
long distances, i.e. without short-distance suppression factors such as 1/Myy. Notice that, 
despite being local, these dimension-eight operators do not appear in the quark effective 
Lagrangian 3 . 

The appearance of these operators has been correctly emphasized in 6 , particularly in 
approaches with the use of a momentum cutoff. However, we remark that this contribution 
from dimension-eight operators is not caused by the cutoff and, in fact, is also true even if 
one uses a purely dimensional regularization, as in Eq. Clearly, a physical effect can 

not depend on the particular regularization employed. What happens in approaches using 
a momentum cutoff A (and does not happen in dimensional regularization) is that, by pure 
dimensional analysis, local dimension-eight operators can also appear at distances > 1/A 
suppressed by 1/A 2 . 4 Moreover, the authors of |Hj have cautioned that these long-distance 
contributions from dimension-eight operators may give potentially large contributions to weak 
matrix elements. 

In our large- -ZV C expansion approach, based on dimensional regularization, we shall see 
i.e. it can be written as a ratio of two polynomials. 

2 If enough conditions were known one could also determine the masses, of course. For a model in which 
this is done, see Ref. U| 

3 In a way, they correspond to having integrated all the meson resonances out of the Lagrangian. This is 
why only Goldstone degrees of freedom appear in JjjJ. This integration can not be encoded in a quark effective 
operator in the Lagrangian because this would give rise to a double counting problem: this quark effective 
operator would generate both Goldstone mesons and the resonances which one thought one had integrated 
out!. 

4 In the lattice jargon they are known as 0(a) effects. For us, a ~ 1/A. 
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that dimension-eight operators produce contributions which can be roughly characterized as 
of 0(a s 8 2 /Hope)i where 5 2 is a Goldstone-to -vacuum matrix element and. i^ope is 

the scale 

controlling the Operator Product Expansion of correlators in pure QCD, i.e. without weak 
interactions. To be a bit more precise about this scale, we may define it in general as the 
momentum scale above which the OPE in pure QCD starts making sense, even if this is only 
in the limited sense of an asymptotic expansion. Phenomenologically this scale /x ope may in 
general be somewhat dependent on the quantum numbers involved, but it is usually of the 
order of a resonance mass (i.e. ~ 1 GeV). 

Besides the long-distance contributions alluded to above, there are other short-distance 
contributions which are also due to local dimension-eight operators but which appear after 
the integration of the W boson or, in general, of any heavy field. These contributions, unlike 
the former, are encoded in the quark effective Lagrangian. Indeed, when writing the effective 
Hamiltonian of Eq. (Q) one is neglecting dimension-eight operators at short distances which 
appear below the charm mass, m c , after integration of this quark. Consequently, these 
dimension-eight operators appear suppressed by the short distance scale 1/m 2 . The size of 
this 1/mj? effect is also an interesting problem in itself but, since it can be dealt with by 
ordinary methods, it will not be the subject of the present work 5 . We plan to come back to 
it in the future. 

In this work we would like to quantitatively investigate the impact of the l//Xo pe -suppressed 
dimension-eight operators in approaches based on dimensional regularization, such as the HA. 
If dimension-eight operators are to give an important contribution in this case, they should 
correct significantly the lowest order term in the OPE at large Q 2 , i.e. the high- z tail of func- 
tions such as PF(z) in the determination of <?as=2(a0 i n Eq. In the following we present 
a new evaluation of Bk using the Hadronic Approximation to large- N c QCD with inclusion 
of dimension-eight operators in the OPE. We now turn to the details of this analysis. 

After chiral corrections, the result in the N c — ► oo limit stems from the full factorization 
of the operator (J3J into two left-handed currents with the result 

B K (N C - oo) = jj and C A s= 2 (m) (N e -»■ oo) = 1 . (12) 

Including the next-to-leading order in the 1/N C expansion one obtains 

3 

Bk = j Cas=2(/^) 5A5=2(/^) , (13) 

where <7as=2(/-0 is given by Eqs. (fTT|) and C&s=2{l i ) is given by Eq. ©. We emphasize 
that the function W LRLR {Q 2 ) in Eq. © is of 0(N C ) and, therefore, it is given by tree-level 
diagrams with infinitely narrow resonances. After the redefinition (|l(Jj) . the function VF(z) is 
of 0(iV c o ) and, recalling that F 2 ~ 0(N C ), #as=2(m) is indeed of the form ~ 1 - 0(l/N c ), as 
one would expect. 

In a combined large--/V c and quark mass {m q ) expansion, the result in Eq. (|13j) is only 
corrected by terms of 0(m q /N c ). Furthermore, it is also interesting to emphasize that if Fq 
is used instead of F% in the definition of Bk in Eq. I|4"]). then the above expression (|13[) is 
also what is being computed on the lattice as the Bk parameter extrapolated to the chiral 
limit. This fact is very convenient as it allows for a comparison between the present large- N c 
method and the lattice results in this limit. 

Figure H shows the different topologies of the resonance diagrams contributing to W(z) 
in Eqs. (|Bllll|) in the large- N c limit. Use of the Mittag-Lefner theorem for meromorphic 
functions 9 allows one to write W(z) as 

5 Some l/m c effects have been afready considered, for exampie, in Ref. 0. 
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Figure 1: Resonance diagrams contributing to W(z). 



^SGi^ + (^ + ot) (14) 

where pi = mf/fj,^^, with the mass of the resonance i; and Ai,B{,Ci are constants. In 
the large- N c limit, the sum is in principle extended to an infinite number of resonances, i.e. 
M — > oo. On the other hand, as mentioned above, the approximation Wha(z) is defined by 
restricting the sum in Eq. (j!4|) to a finite N . In particular, the work of Ref. [2j restricted 
J\f to just one vector resonance, whose residues where matched onto the leading and next-to- 
leading terms in the chiral expansion but only to the leading (nontrivial) term in the operator 
product expansion of W{z). This way it was obtained that Bk = 0.38 ± 0.11; with the error 
being an estimate of uncalculated terms, including 1/N C corrections. The purpose of the 
present work is to compute the next-to-leading term in the operator product expansion of 
at large-z, and its impact on the previous value for Bk- 

2 Dimension 8 operators in the weak OPE 

Let us consider the Green's function ©. The large-Q 2 behavior of this function is governed 
by the Operator Product Expansion 

/ otx e «*» r{L*(*)Z&(0)} = E c S 6) (Q 2 )0? } +£ c ? ) (Q 2 )Of ) + , (15) 
J i i 

where the first (second) sum runs over a set of dimension six (eight) local operators to be 
discussed in the following. In the case of dimension-six operators this set is actually limited 
to just one, to wit 

0(6) = s L ^d L (0)s Lltl d L (0) , (16) 
with the Wilson coefficient given, to lowest order in a s = gl/^, by 

c(6)(Q2) = .i^ (17) 

The relevant diagrams are listed in Fig. [2j 

The determination of the dimension-eight operators is more involved. We used the 
Schwinger's operator method and obtained (after use of the equations of motion) 

0f ] = s^^O(0) sT u a d(0) + STtPpWdfi) sT u a d(0) + 
sT»d(0) I^^rXO) + sT v a d{0) sT a u V^d{Q) , 
Of = sO^(O) sT a v V^d(0) + s^r-d(O) sV~^ a d(0) , 
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Figure 2: Diagrams contributing to the dimension 6 terms in the OPE in Eq. 115\) . 



of = s%Kd(0) sT a u V^d(0) + sT a u V^d(0) 8^1^(0) , 

of = s%vy(Q) s^r^(o) + sr a u v^d(o) sT^d(o) , 
of = i^r^(o)ir^d(o) + ir^(o)i^r^(o), 

of = g s G%(0) {sT%d(0) sT»d(0) - ST"d(0) ST£d(0)} (18) 
where the following conventions were adopted: 

Gfj,v = d^-d u A^-ig s [A^A„] , G^ = ^G^e^ , with e 0123 = +1 . (19) 
The corresponding Wilson coefficients are 

<f = ^> f f, (20) 



where 



(8) 5 ( 8 ) 22 ( 8 ) 8 

m = 3 . % = y i ^ 3 

r7 (8) - — r? (8) - — r7 (8) - — (21) 

This result agrees with the calculation in Ref. [B] if one takes into account that, in the present 
case, the specific flavor structure (sd) 2 leads to some simplifications. In particular, it allows 
one to rewrite Q4 in Ref. H in terms of the operators in Eq. ((TH|) . 

Inserting the operators (|16I18|) in Eqs. (|15l6llllj) one can compute the expansion in powers 
of l/z of the function VF(^) in Eq. (|llj) . The terms of 0(z^ 1 ) were already computed in Ref. 
[2]. In doing the calculation of the next-to-leading terms, of 0(z~ 2 ), we notice that several 
further simplifications take place. Firstly, since the contribution is proportional to a s - see 
Eq. ()17I2U() - , the large-iV c limit allows one to use factorization. Secondly, the contribution 
from Oq is 1/N C suppressed due to the value of its Wilson coefficient, Eq. (|2*Tj) . Thirdly, one 
is only interested in the leading term in the I — > limit of the Green's function Q. Since, 
after Fierzing, one ends up with matrix elements of the form 



< 0|s L r^ M d L (0)|^(/) >~ag^ + b l^ly , (22) 

by contracting with g^ v one immediately concludes that a ~ 0(l 2 ) and, consequently, 02,3,4,5 
can be neglected since they yield contributions to Eq. © with an extra power of I 2 . Therefore, 
the term of 0{z~ 2 ) in the large-z expansion of W(z) in Eq. (jll|) is governed solely by 
0\. Furthermore, using D 2 = ft 2 + ^-cr^G^ in the expression for Oi, and the previous 
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considerations, lead to the conclusion that the final result can be given in terms of a single 
matrix element 6 , namely 

< 0\g s s L G» v \ a ^d L \K(l) >= -iV2F b\ i„ , (23) 

where b\ is a parameter to be determined in the next section. In conclusion, the OPE of 
W(z) in Eq. (flip can be written in terms of this parameter 5\ as 

Since the term proportional to b\ yields an ultraviolet convergent integral in Eq. pi|). it can 
be computed in 4 dimensions. This is why there is no e dependence in Eq. (]24[) accompanying 
this term. 



wOPE{z) = 2^a s Fj 

Vhad 



3 Determination of b\ 

The mixed quark-gluon matrix element (|23j) was determined in the literature in j 1 01 lllj . 
Both approaches relied on an analysis of Borel Sum Rules applied to two-point functions 
with nonvanishing contribution from the perturbative continuum. Although their final result 
for b\ happened to be in agreement, there were some criticisms raised by on the analysis 
of ^0], in particular on the choice for the onset of this perturbative continuum, sq. This 
partly motivated us to do a reanalysis of the matrix element Q23JI . 

A particularly simple way to bypass the need for choosing a value for sq consists in looking 
at Green's functions which are order parameters of spontaneous chiral symmetry breaking 
since in these functions the perturbative continuum vanishes by construction. This procedure 
has proved fruitful in several contexts ^2] and will be applied to the present case as well. 

Consider the Green's function 

n^(Q 2 ) =ij d 4 x e iqx < \T^s L G a ^ a d L (x) d Rlu s R (0)^ |0 > . (25) 
Lorentz invariance guarantees that, in the chiral limit, it has the following tensorial structure 

^(q) = (q,qu-9,uq 2 )fl LR (q 2 ) ■ (26) 
The large-g 2 fall-off of the function (|26|) can be straightforwardly computed to be 

n l V) = -f^|^ , Q^-i, (27) 

where factorization of the four-quark condensate in the large- N c limit has been used. Since 
N-lr obeys the unsubtracted dispersion relation 

CO 



n LR m = I dt~ -;"^" J , (28) 



HmU LR (t) 
'tt t + Q 

we can consider as the minimal hadronic approximation the following spectral function 



o 



-ImU LR (t) = -^5(t) + " ml) (29) 

TT 4 8 



6 This is related to the fact that sj^G^d is the only nontrivial dimension-five operator in the chiral limit 
with these quantum numbers |f()| . 
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Figure 3: Plot of the right and lefthand side of Eq. \35}) (in arbitrary units), corresponding to the 
dashed and solid curve, respectively. 



in which one introduces, besides the Goldstone boson, a vector resonance. Inserting 1)28)) into 
and expanding for large Q 2 one obtains, upon comparison with the OPE of (J27J), the 
two Weinberg-like sum rules 

4 8 

= f as <^>\ (30) 

from which the unknown 5 V and b\ parameters are readily determined to be 

2 16ir a s < > 2 



5 



V 9 f v m\ 



Using F ~ 0.087 GeV, Jy ^ 0.15, my ~ 0.77 GeV, a s {2 GeV) ~ 1/3 and (^)(2 GeV) = 
-(280 ± 30 MeV) 3 one obtains 7 

5% = 0.12 ± 0.07 GeV 2 , S v = 0.06 ± 0.04 GeV 4 . (32) 

Strictly speaking, both parameters 6^ and 5 V depend on the renormalization scale fi. This 
dependence, however, is very small. As a matter of fact, according to Ref. [T3] the combina- 
tion 

a s (^y~ QsSLl^G^dL (33) 

is renormalization group invariant. This means that the [i dependence of the parameters 
5k y in Eq. 1)3 ljl yields a variation of 10% if \i is varied in the range 1 < fi < 2 GeV. 

Apart from the dependence on the quark condensate, in order to check how much the 
result (|32|) depends on the two-particle decomposition of the spectral function in Eq. l)29]) . 
we can now repeat the analysis by introducing a third state. We take this state to be the 
first axial resonance. In this case, the spectral function reads 

-Imfl LR (t) = 5(t) + ^5(t - m v ) + %t - m\) . (34) 
7r 4 8 8 



7 We have added generous error bars in the quark condensate to include the present spread of values in this 
quantity U^J. This error is the dominant one. 
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One can now apply the Borel transform to (|34[) and the OPE of (|27|) to get the relation 



_ ^pL + ZM e -^r + fv^^r = < ^ >2 T ^ (35) 

4 8 8 9 

where r is the Borel parameter. The equality above is fulfilled within a rather wide window 
of duality, namely < r < 1, as can be seen in figure Q- Using /a = 0.08 and tha = 1.2 
GeV, from this plot one extracts values for b\ , S v which agree with the ones obtained in 
(j32|) . There is, therefore, consistency in the result obtained with the two methods. For the 
new parameter, S\, one obtains ~ 0.05 GeV 4 . We remark that the value of 5\ is a bit lower 

than the results obtained in refs. |10| lllj. To the best of our knowledge there is no lattice 

•2 

A- 



evaluation of S 2 



4 Numerical Evaluation of Bk and Conclusions 



In Ref. [2] the correction from dimension-eight operators was not considered and matching 
to the OPE was accomplished with just a single vector resonance whose mass was taken to be 
the physical p mass. In the present case we have a further constraint given by the new large- z 
behavior in Eq. (|24|) caused by the presence of the parameter b\. Although one might think 
that this new constraint could allow one to determine the mass of the vector resonance, it 
turns out that this is actually not possible: the positive value of 5\ would force the p mass 
to be purely imaginary. Therefore one must consider a further resonance. 

Besides vectors, scalar particles also contribute to the coupling £3^3^!]. The dependence 
of the function W(z) on L3 at low-z found in |2j -see Eq. (|37jl - suggests, therefore, a scalar 
particle as the new resonance in the sum (|14|). Furthermore, a scalar particle appears as a 
single pole in the function W(z) in Eq. (|14|) as a consequence of the quantum numbers being 
exchanged in the Green's function © (see Fig. 1). Interestingly, one scalar residue is all one 
needs to balance the new constraint from the OPE encoded in b\ through Eq. (|24j) . As to the 
scalar mass, which remains undetermined, we allow for the generous variation m s = 900 ±400 
MeV. 

Gathering all the pieces, one has as the interpolating function 



W HA (z;{S,V}) 



Ay 



By 



+ 



+ 



As 



(z + pv) (z + p v ) 2 (z + py) 3 (z + ps) 



(36) 



to be matched onto 

W* PT (z) 



2 

6 _ 24 ^f* (2Li + 5L 2 + L 3 + L 9 ) z + O (; 



F 2 
^0 



(37) 



at low energies, and 



W 



OPE, 



z) 



24na s F^ 1 

Vhad Z 



si 



+ 



(38) 



at high energies. This results in the following 4 constraints: 



Ay + As 



By - A v pv - Asps 



247ra,F 2 

^had 

247ra s F 2 

^had 



1 + ^(5 + ^) 



10 2 
9 6k 
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Figure 4'- Plot of the large- and small-z expansions of the function W(z) in Eq. as given by the 

OPE and Chiral Theory, respectively (dashed curves). The solid curve corresponds to the interpolating 
function Wha(z) obtained in Eq. |5*61) . 



Ay By Cy A S 

1 S- H n 1 = 6 

Pv Py Py PS 

4^ + 2^ + 3% + ^ = 24%(2L 1 + 5L 2 + L 3 + L 9 ) , (39) 
Pv Pv Pv Ps F o 

for the 4 unknowns As, Ay, By and Cy. The combination (2L\ + 5L 2 + L3 + -^9) in the last 
equation equals 11.2 x 10~ 3 if its experimental value evaluated at the p mass scale is used 
|17j 8 . This combination is renormalization scale dependent, but this dependence can be 
neglected as it is subleading in the 1/N C expansion. Numerically, it amounts to 30% if is 
varied in the range 0.5 < fi < 1 GeV. The coupling a s is evaluated at the scale Hhad, which 
is chosen in the range 1 — 2 GeV. All these variations will be included in the final error. 

Figure 4 shows the low- and high-z behavior of the function VF(z) given by its chiral and 
operator product expansion (dashed curves) together with the interpolator Wha(z) in Eq. 
(f3*oT) (solid curve). 

The main effect from the dimension-eight operators is to push the large-z tail of the 
function W(z) slightly upwards, as a consequence of b\ being positive. This implies a larger 
value of the integral in Eq. (|lljl . i.e. a larger area under the curve and, consequently, a 
smaller value for Bk- One can get a rough idea about the size of the contribution from the 
dimension-eight operators by noticing that the 5 2 K in Eq. (|38j) generates a correction to the 
coupling gAS=2 in Eq. (|TT|) given by 

l4 ad 24vra s (^^)F 2 10 b\ f°° dz 

ogAS=2 ~ 



(40) 



had 



18 vr nl 



"ope 

where /j, ope is the scale above which the large-z expansion (i.e. the OPE) starts making sense. 
8 If the results of Ref. are used, this combination is 8.7 x 10~ 3 (11.7 x 10~ 3 ), respectively. 
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Numerically, one obtains \5g\ < 0.03 when fi ope ~ 1 GeV and the value for 8 2 K in Eq. (|32j) 
are used. 

To get a more accurate result one has to use the solution for As, Ay, By and Cy from 
Eqs. (|39|) to construct the hadronic approximation Wha m Eq. (|36|) . One can then explicitly 
perform the integral in Qll|) to obtain, after multiplication by the Wilson coefficient in Eq. 
(J2J) according to Eq. a cancelation of scale and scheme dependence with the following 
result 9 



Using now the values a s (2 GeV) ~ 0.33, F = 0.087 GeV, 1 GeV < fi had < 2 GeV, 
the masses my = 0.77 ± 0.03 GeV, ms = 0.9 ± 0.4 GeV and the matrix element 5 K = 
0.12 ± 0.07 GeV 2 , the expression (J3TJ) yields 



where the error quoted embraces an estimate of uncalculated 1/N C corrections and the "noise" 
created by the uncertainty in all the parameters, taken one at a time. 

We would like to conclude with a comparison with other results. Since the AS = 2 
operator © sits in the (27l,1r) representation of flavor SU(3)l x SU(3)r, the result (|42|) 
translates into a value for the coupling constant 527 which governs A/ = 3/2 processes, such 
as — ► 7T + 7r°, namelv |2()l I21j 



This number is in very good agreement with a recent extraction [22] of 527 from K — » 37r 
decays which yields, after subtraction of chiral corrections (within reasonable assumptions), 
the value 527 — 0.24. 

However, because chiral symmetry is much more difficult to have on the lattice than on 
the continuum ,23] , the situation concerning the value of Bk in numerical simulations in the 
chiral limit is not totally clear. It is seen in lattice data that Bk dips as quark masses go 
to zero|24j. but the errors do not yet allow a sufficiently accurate extraction 10 . Nevertheless, 
we remark that two lattice collaborations have recently obtained Bk — 0.3 — 0.4 when 
extrapolated to the chiral limit 25 . 

Finally, alhough Bk in the chiral limit may be an interesting problem per se, it is clear that 
Nature is not in this limit and, therefore, that it is of the utmost importance to compute chiral 
corrections to the result (|42|) . As a matter of fact, lattice results |27j and phenomenological 
analysis [2B] seem to favor a large value for Bk at the physical kaon mass of the order of 
twice the value in (|42j) . This would imply that the 0(m q /N c ) corrections to our result ()42j) 
are very large; a result which should be understood by analytic methods. That this is not 
impossible was shown in the calculation of Ref. [2H]. Even though within a model with 
notable differences with respect to QCD |30| . the authors of Ref. [22] obtained a value of 
Bk in the chiral limit in agreement with our result (|42[) while at the same time producing a 
larger Bk at the physical kaon mass, in agreement with Refs. [271 128j . 

9 The quadratic dependence on Hhad is actually fictitious as is canceled by that of Av,s,Bv and CV; see Eq. 
10 See Ref. for a new strategy to try to overcome the difficulties. 




(41) 



B K = 0.36 ±0.15 



(42) 



g 27 = -B K = 0.29 ±0.12 . 
5 



(43) 
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